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Quark Coalescence and Charm(onium) in QGP 
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Abstract. The potential of heavy quarks as probes of the environment produced in hadronic and heavy-ion 
reactions is discussed. A key role is played by coalescence processes and/or resonance formation which are 
promising candidates to provide a comprehensive understanding of phenomena associated with reinterac- 
tions of both open and hidden heavy-quark states. 

PACS. 12.38.Mh Quark-Gluon Plasma - 25.75.-q Relativistic Heavy-Ion Collisions - 14.40.Lb Charmed 
Mesons 



1 Introduction 

In hadronic and heavy-ion collisions, heavy quarks (Q=c, 
b) are believed to be (almost) exclusively pair-produced 
(QQ) upon first impact in hard partonic collisions yQ. 
This renders them excellent agents of the subsequently 
formed medium and their reinteractions within. The lat- 
ter include: (a) coalescence with surrounding quarks as 
hadronization mechanism in addition to fragmentation, 
thereby probing the chemical and kinematic properties of 
the medium [2 3,4 5,6]; (b) energy loss of high- momentum 
Q-quarks 0IE]> which, with increasing interaction strength 
toward lower momentum, eventually leads to (c) thermal- 
ization [9*lll0lllllll2| : and, if the latter can be established, 
(d) in-medium dynamics of open and hidden heavy-flavor 
states 13,14,15.16,17 , which is particularly exciting in 
view of recent QCD lattice calculations 18,19,20 indi- 
cating the survival of low-lying charmonia well into the 
Quark-Gluon Plasma (QGP). In this paper we will ad- 
dress the above issues essentially in that order. 

2 Coalescence in Hadronic Collisions 

In elementary hadronic reactions (pN, ttN) evidence for 
reinteractions of c-quarks arises from (large) flavor asym- 
metries in D-meson production yields. The asymmetries 
are most pronounced at forward rapidities (or xp), suc- 
cessfully being attributed to coalescence of c-quarks with 
valence quarks of the projectile 4 . The pertinent recom- 
bination cross section can be written as [2] 



dxp 



diXq 



dz 



d 2 a (c ^ 
dxgdz 



1Z(x q , z; xf) 



where the main elements are: (i) the c-q production cross 
section composed of a 2-parton distribution function (2- 

(2) 

PDF), (where i=g,q,q participates in the hard process 



to produce the cc pair), and the standard perturbative 
QCD (pQCD) cc cross section, and (ii) the c-q-^D recom- 
bination function, 1Z. The 2-PDF is usually factorized into 
two single PDFs with phase space correction, 



fg> = cf q (x q ) Mx t ) (i-x q - Xl y , 



(2) 



whereas 1Z represents a Z?-meson wave function which in 
Ref. [H] has been assumed to be Gaussian in rapidity space, 



ft(2/,-,2/c,y) = exp(Zh/ 2 /2^)/ 



(3) 



This form of the recombination function [22 allows to gen- 
eralize the coalescence formalism to include sea-quarks 6 , 
and thus address flavor asymmetries also at central xp, cf. 
Fig. H The experimentally observed asymmetries in inclu- 
sive yields (xf>0) are quite appreciable, e.g., D~ /D + = 
1.35±0.05 (versus 1 in isospin-symmetric fragmentation), 
D°/L»°=0.93±0.03 (vs. 1) and L> ± /( J D°+D°)=0.415±0.01 
(vs. 0.33) for fixed-target n~ N collisions (averaged over 
a weak energy dependence for v / s=19-34GeV) [231 - T ne 
data are rather well reproduced by a combined coalescence 
+ fragmentation approach [5j (for a somewhat different 
framework based on power corrections, see Ref. [H])- 



3 Open Charm in the QGP 

Final-state interactions of heavy quarks arc enhanced when 
embedding them into a heavy-ion collision, where, at ul- 
trarelativistic energies, intense reinteractions of light par- 
tons are believed to form locally thermalized matter within 
a time of r<lfm/c. At high momenta c-quarks rescatter 
perturbatively inducing a softening of the primordially 
power-like pt-spectra, with subsequent hadronization in 
the vacuum (fragmentation). The predicted suppression 
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Fig. 1. D-meson flavor asymmetries, A=(JVo 1 — iVo 2 ) /(Nr> 1 + 
Nd 2 ), in n~-A reactions |22I23) compared to coalescence + 
fragmentation calculations Upper data points and curves 
are for Di~D~ and D2—D + , whereas lower data points 
and curves are for Di=D° and D2=D°. "Leading" particles 
(D~—cd, D°=cu) are defined as sharing a valence quark with 
the projectile (n~=du); note that the approximate absence of 
an asymmetry for D° /D° is accounted for in the model due to 
the predominant production of cc pairs in the forward direc- 
tion via uu annihilation, rendering the valence u unavailable 
for recombination. 



factors relative to p-p collisions range from 0.2 [H] to 0.5 0, 
with a rather small azimuthal asymmetry, V2<5% |24) . 

Toward lower p t , the phase space density of the medium 
increases and coalescence with light quarks is expected 
to become competitive 25 26 2Z] - The same expression, 
Eq. Q), can be applied with the light quark distributions 
being replaced by thermal (+ quenched pQCD) ones as es- 
tablished from light hadron production systematics (also, 
the recombination function 1Z is typically substituted with 
a hadron wave function in transverse momentum). The 
extension to low p t is, in principle, more controlled than 
for light-light (q-q) coalescence, since at the scale of the 
hadronization temperature, secondary c-production is neg- 
ligible. With previously determined light-quark distribu- 
tions, charmed-hadron spectra become a sensitive probe of 
the dynamics of c-quarks in the QGP. This has first been 
quantified in the context of "charm-like" single-electron 
spectra in Ref. [2HJ, showing that u|(Pt) : (a) closely reflects 
the V2 of the parent D-meson, (b) exhibits a marked dif- 
ference of more than a factor of 2 between the cases where 
the c-quark distributions are either taken from p-p colli- 
sions, or assumed to follow the systematics of light quarks 
(including collective expansion) , cf . Fig. [21 Current data 
at RHIC from PHENIX [2H| and STAR (preliminary) (2H| 
seem to favor the quasi-thermalized scenario. If confirmed, 
this raises at least two further questions: 
(i) Is the predicted i>f consistent with pertinent p t -spcctra 
(i.e., the ratio of central Au-Au to collision-scaled p-p 
spectra, Raa)1 
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Fig. 2. Azimuthal asymmetry of "non-photonic" single-e 
spectra in minimum bias Au-A u(y/s jvjv=200GeV) |28I29| com- 
pared to coalescence model predictions |25| using c-quark dis- 
tributions from either p-p collisions (t>2=0, solid line) or as- 
suming a transverse flow and «2-profile as determined for light 
quarks from fits to light-hadron spectra (i;2>0, dashed line). 
The band indicates predictions from jet-quenching |24| appli- 
cable at sufficiently high p t . 



(ii) What are microscopic mechanisms for thermalization 
of c-quarks (or D-mesons)? 

Concerning (i), it has been pointed out [23] that for 
single-e 1 * 1 pt-spectra in central Au-Au, 7V co /;-scaled D-me- 
son spectra from p-p collisions lead to results rather similar 
to a scenario based on full thermalization and collective 
flow close to hadronic freezeout (T~130MeV, v^=0.65), 
due to large blue shifts with mD = 1.87GeV (also, bottom- 
decay contributions become significant above pj~3GeV). 
However, in hydrodynamic analysis coupled with a 
Fokker-Planck treatment of c-quarks in the QGP, a v 1 ^ 
of -10-15% is associated with R AA (p t >3GeV)<0.1. Coa- 
lescence model calculations [25], based on recombination 
at the phase boundary, imply a suppression factor sim- 
ilar to jet quenching, R AA (p t ~3GeV)— 0.2-0.5, but with 
t>.P(p t ~3GeV)~15%. Note that when starting from c-quark 
spectra, fragmentation leads to a degradation, whereas co- 
alescence to an increase, of the resulting D-meson p t . 

Concerning (ii), it has been known for a while JjJ] (and 
confirmed in Refs. |10llllp that perturbative c-quark rescat- 
tering off quarks and gluons in the QGP implies kinetic 
relaxation times r c " lerm >10fm/c for T~400MeV, too long 
to achieve thermalization at RHIC. However, as shown re- 
cently nonperturbative rescattering in the QGP can 
lead to a substantial acceleration of equilibration: imple- 
menting the notion of D-meson-like resonances within a 
Fokker-Planck equation, a reduction of r *' lerm by a factor 
of —3 as compared to using pQCD cross sections has been 
found (for T<2T C ), cf. upper panel of Fig. [3] (similar for 
6-quarks, but with absolute values T t herm ^T* herrn ). The 
main difference in the two mechanisms resides not so much 
in the total cross sections (lower panel of Fig.[3J), but in the 
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Fig. 3. Upper panel: c-quark equilibration times in QGP with 
pQCD interactions (upper band; upper [lower] line corresponds 
to a s =0.3[0.5])) and when adding "D"-meson resonance rescat- 
tering (lower band; the range of resonance widths indicates 
variations in the coupling constant of the c-q-D vertex with 
the upper [lower] line corresponding to Ji5=0.3[0.8]GeV) 
Lower panel: underlying total c-parton cross sections. 



isotropic angular distribution for the resonance case as op- 
posed to forward-dominated pQCD scattering. It has also 
been noted that the efficiency of this mechanism relies 
to a significant part on the D-states being located above 
the c-q threshold (i.e., not being boundstates, which ren- 
ders them inaccessible in 2^2 scattering, especially due to 
the thermal energies carried by the light quarks) . It will be 
very valuable to check this in QCD lattice calculations, as 
well as whether previously found q-q and Q-Q states carry 
over to the Q-q sector. Furthermore, an increasing popula- 
tion of (colorless) "hadronic" states in the cooling process 
toward T c could serve as a mechanism to put phenomeno- 
logically successful coalescence models on a firmer basis 
(also in the light-quark sector). The in- medium mass of 
open-charm states in the QGP also bears on the produc- 
tion of charmonia, as will be seen in the following section. 



From a phcnomcnological point of view, it should be 
kept in mind that any process contributing to elastic c- 
quark scattering in the QGP, c + X\ — > c + X%, in prin- 
ciple also gives rise to secondary cc production in the 
crossed channel, X\ + Xi — > c + c, which can be con- 
strained experimentally by total cc yields (including non- 
trivial centrality dependencies). E.g., in Ref. [27J it has 
been pointed out that when upscaling the perturbative 
gc — ► gc cross section by a factor of 3 (to generate an el- 
liptic flow comparable to light quarks), secondary charm 
production is at the 40-50% level of the primordial yield 
in central Au-Au(y/s=200AGeV). This is expected to be 
less pronounced for heavier exchange particles, such as 
"D" -mesons. 



4 Charmonium in the QGP 

A central quantity in evaluating medium effects on quarko- 
nium states, \P, in a heavy-ion collision are their inelastic 
cross sections, a^ lss , with partons in the QGP, determin- 
ing the pertinent dissociation rate as 

d 3 k 



(Wfe,T) V re l (T^ SS (s) . (4) 



A widely used model for er* ss is the gluon-absorption 
break-up |3"Dl3"2*] , g + 9 — > c + c, characterized by a pro- 
nounced maximum at a gluon energy u> max ~l.5eB (f-B- 
quarkonium binding energy), see lower panel in Fig.0J For 
J/tp mesons with their free binding energy, e B ac =6A0MeV, 
u m ax essentially coincides with thermal gluon energies, 
lu='3T, for T~300MeV. Debye screening of the Q-Q poten- 
tial in the QGP is, however, expected to substantially re- 
duce cb |34| . This renders gluodissociation an increasingly 
inefficient process at higher T due to a shrinking break- 
up kinematics, cf. dotted lines in Fig. For small es, 
"quasifree" dissociation [55I1T3] . a(a. q)+ x I / — > c+c+g(q, q), 
albeit formally suppressed by one power of ce s , has been 
identified as a more important mechanism due to much 
larger overlap with the thermal (quark + gluon) phase 
space (cf. solid and dash-dotted lines in the lower panel of 
Fig. ID. 

If the number of heavy quarks in a heavy-ion collision 
is large enough, their recombination into quarkonia could 
become a significant (or even dominant) contribution to 
the final yield 35 36 37 33 38 39 . The conditions for this 
to happen can be assessed in terms of a simple rate equa- 
tion for the time evolution of the number of tf^'s, 



dt 



= -r 9 (N* - N e w q ) 



(5) 



Besides the reaction rate TV, the other quantity govern- 
ing the evolution of Nq, is the equilibrium abundance, 
Np(T; 7 C ), which determines & regeneration, i.e., the gain 
term in Eq. (J5J, as required by detailed balance. N^(T; j c ) 
is typically evaluated in the canonical ensemble with the 
total number of (primordial) cc pairs fixed via a fugac- 
ity 7c=7e=e Mc ^ T . This implies that N^ q (T; j c ) is sensitive 
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Fig. 4. Upper panel: J/ip lifetimes in the QGP using glu- 
odissociation |32| with vacuum (dashed line) and in-medium 
reduced (dotted line) binding energy, as well as quasifree dis- 
sociation 1331 with in-medium reduced binding energy (solid 
line). Lower panel: pertinent cross sections (line identification 
as in upper panel) relative to thermal parton distribution func- 
tions (dash-dotted line). 



to the open-charm spectrum, in particular (in-medium) 
masses of c-quarks (or D- mesons) |15U16j : e.g., if m* (or 
m* D ) is reduced in matter (with constant), c- and 
c-quarks are thermally favored to occur in open-charm 
states, thus reducing N^(T\^ C ). Finally, the gain term 
depends on the c-quark momentum distributions; its par- 
ticularly simple form in Eq. based on thermalized o 
quarks, illustrates the impact of c-quark rescattering (as 
discussed in the previous section) on charmonia. There- 
fore, thermalization of c-quarks opens the window on equi- 
librium properties of both open and hidden charm, i.e., 
their masses encoded in 7V^ 9 (T;7 C ), as well as charmo- 
nium widths A sensitive observable to distinguish 

direct and regenerated J ftp's turns out to be their elliptic 
flow, v% )40|I41I125II42| . If only suppression is operative, wf 
reaches a maximal value of ~2-3% while it grows up 
to ~15% at pf ~4GeV for thermal c-c coalescence |25| . 
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Fig. 5. Time evolution of the J ftp abundance in central Au- 
Au(^/SNN~'200GeV) based on a solution of the rate Eq. JjsJ in 
an expanding thermal fireball including in-medium effects on 
both open and hidden charm mesons |lt)| . Dashed line: sup- 
pression only (no gain term); dash-dotted line: temperature- 
dependent equilibrium number, NjJ^(T;*/ c ); solid line: total 
number Nj/^(t). 
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Fig. 6. Theoretical predictions for the excitation function 
of J ftp production in central Au-Au collisions. Suppression 
(dashed line) and regeneration (dash-dotted line) components 
combine into a rather flat energy dependence for the total yield 
(solid line) E3- 



A calculation ^H] of the time evolution of Nj/^ in cen- 
tral j4it-Au( v / s=200AGeV) based on Eq. JSJ) including in- 
medium masses of open-charm and reduced J/\P binding 
energies, as well as incomplete thermalization of c-quarks 
in the early stages, is displayed in Fig. 03 One finds that 
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Fig. 7. Theoretical predictions for the centrality dependence of 
J/ip production in intermediate-size-ion collisions at SPS (up- 
per panel) and RHIC (lower panel) , including both suppression 
and regeneration processes |16) . The bands in the upper panel 
reflect uncertainties in the formation time which is expected 
to increase at lower collision centrality (dashed lines). 



the J/ip yield equilibrates close to the phase boundary, 
with the major contribution arising from regeneration in 
the QGP and little changes in the "mixed" and hadronic 
phase. Note that this result crucially hinges on the notion 
of the J/tp surviving as a resonance in the QGP under 
RHIC conditions, T<2T C . The final yield is a factor of ~4- 
5 increased over a scenario with suppression only. The sit- 
uation is quite different at SPS energies (y/s—17 .3AGeV): 
with an expected open-charm number iV cg ~0.2 in central 
Pb-Pb, secondary charmonium formation is negligible and 
J/ip suppression is the main mechanism at work. Obvi- 
ously, this calls for mapping out the excitation function 
for y / s=20-200GeV (accessible at RHIC), as suggested in 
Ref. 33 ■ Based on Fig. one expects a transition from 
a suppression-dominated regime (SPS or low RHIC ener- 
gies) to a regeneration-dominated one at -y/sXLOOAGeV, 



resulting in a rather flat energy dependence (possibly with 
a shallow minimum). 

Complementary information on the interplay between 
primordial and secondary J/ip production can be extracted 
by going to smaller nuclear collision systems. Pertinent 
predictions are shown in Fig. reconfirming the absence 
of noticeable regeneration at SPS (as well as a smooth cen- 
trality dependence; upper panel), but an approximately 
equal amount of primordial and regenerated J/ip's for cen- 
tral Cu-Cu at RHIC (lower panel). 



5 Conclusions 

Hadrons containing heavy quarks are excellent probes of 
the environment formed in nuclear reactions. Evidence for 
coalescence mechanisms in elementary hadronic reactions 
finds its natural extension for both D-mesons and char- 
monia to heavy-ion collisions. In addition, at RHIC, the 
produced medium appears to interact strongly enough to 
thermalize c-quarks (but not 6-quarks). If confirmed, "£)"- 
meson resonance formation in the QGP (coupled with per- 
tinent coalescence at T c ) might be the key to a simultane- 
ous understanding of (suppressed) pt spectra and (large) 
elliptic flow of D-mesons (and single electrons). The tran- 
sition into a perturbative energy-loss picture could be shif- 
ted to higher pt than for light hadrons. Resonance states 
in the QGP also have substantial impact on charmonium 
production, facilitating their regeneration in the 1-2T C re- 
gime where inelastic collision rates are high. Here, ther- 
malization of c-quarks would enable a rather direct win- 
dow on spectral properties of open and hidden charm, i.e., 
their masses and widths. Work in progress on T produc- 
tion 03j seems to indicate, however, that even at LHC 
their suppression is prevalent, due to a lack of thermaliza- 
tion of bottom quarks. Thus, a simultaneous observation 
of T suppression and the absence thereof for J / ip at col- 
lider energies would provide strong evidence for secondary 
charmonium production. 

Among the main challenges yet to be met is establish- 
ing connections of heavy-quark observables to ( "pseudo"-) 
order parameters of the QCD phase transition 01]. With 
low- lying charmonia possibly surviving up to 2T C , their 
dissolution evades a direct relation to T c . A suitable quan- 
tity could be their inelastic width, which in model calcula- 
tions is typically quite different (smaller) in the hadronic 
compared to the QGP phase JH|- Quenched lattice cal- 
culations |2Uj indicate a similar trend, but unquenching 
has to be awaited for more definite conclusions. We also 
mention the recent work of Ref. 45 j , where an increase in 
transverse-momentum fluctuations of open-charm states 
has been linked to a first-order transition. 

Looking into the future, it seems that the combined 
experimental and theoretical analysis of heavy-quark ob- 
servables in ultrarelativistic heavy-ion collisions is on a 
promising path toward providing a milestone in the iden- 
tification of the QGP. 
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